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AIRPLANE STABILITY CALCULATIONS , 
and Their Verification by Flight Tests.* 

By Augusto Rota, 
Italian Naval Engineer. 

Questions relating to the stability of airplanes in flight 
have greatly interested theorists since the beginning of aviation 
and there is no one who is not acquainted with the remarkable 
works of Mr. Pairileve, Prof. Bryan, Col. Crocco and many, others. 

Unfortunately nearly all these works were only veryremarka- 
ble solutions of a very interesting theoretical mechanical problem, 
but of very little importance from the practical point of view. 

In fact, constructors proportion their airplanes, in their 
organs of equilibrium and stabilization, according to practical 
data acquired by years of observation and successive improvements 
of old-.;, airplanes, with little thought of the results of the dis- 
cussion, of the systems of differential equations established by 
scientists. 

For some time, however, they have begun to occupy themselves 

with questions connected with longitudinal stability. As was qv.it e 

natural, they have sought to divest the theories of everything 

complicated, in order to employ the simplest and also the surest 

data. Thus, in a manner similar to that followed by the engineers 

of naval construction, they have attached the greatest importance 

to the coefficient of initial longitudinal stability, i.e., to the 

quantity ^ ^ 

°s - 5 a 

♦From "Premier Congres International de la Navigation Aerienne," 
Paris, November, 1921, Vol. I, pp. 42-48. 
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in which they consider 

M = M (V a 6, a, t, qp f ). 
M being the moment with respect to the center of gravity of the 
airplane at the speed V, the relative density 6 the angle of 
attack a the angle of deviation of the elevators t and the an- 
gle " of flight <p . • , 

It was very natural to consider also the Quantity 

H -. tt 6a 

tt being the weight of the airplane, H the "virtual metacentric, 
height" and 



h = 



H 



V 3 5 

h receiving the designation of "reduced metacentric height. " 

The terms 0 , H and h must have positive vaD.ues, so that 
the airplane will be stable in its longitudinal motion. The con- 
sideration of the term H is justified by analogy with the cor- 
responding term in the theory of the ship, that of the term h 
in the fact that the. action of all the organs of control and equi- 
librium of the airplane is proportional to the term V 3 6. 

To all these quantities however there has hitherto been gen- 
erally given only a qualitative importance,* not only because it 
was not yet very well known what value to give them in comparison 
with the other characteristics of an airplane, but also because 



* We know, for example, the. diagrams employed by Mr. Eiffel and the 
displacements of the center of gravity for different values of the 
angle t. These diagrams answer quite well for the qualitative 
solution of the problem of stability. 
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the deduction of the results oi 1 experiments on models of the val- 
ues of 

6M 
6 a 

is not made very easily, especially when the v.ceition of the cen- 
ter of gravity has not yet been fixed. 

It was for the purpose of eliminating this second disadvant-.,: 
age that we thought of utilizing, for the quantitative study of 
the stability of airplanes, a diagram which is very simple to con- 
struct ffrom the data of the tunnel tests and which neither de- 
pends on the position of the center of gravity nor of the angle 
of deflection of the elevators. This diagram is constructed by 
means of straight lines drawn through the metacenters of the com- 
plete airplane, in a direction parallel to the tangents to the 
polar of the airplane relative to a system of axes fixed with ref- 
erence to the airplane (Fig. l). 

If m is the metacentric curve of the airplane correspond- 
ing to a given fixed position of the elevator, if p is the po- 
lar with reference to a system of axes fixed with respect to the 
airplane, and if M is the metacenter corresponding to an angle 
of attack a, to which the slope of flight <p , the total action 
of the air on the airplane F is on the straight line tangent at 
M. to the metacentric curve, parallel to the corresponding vector 
OK of the polar. 

If the angle a, without changing the speed nor the direc- 
tion of flight, becomes a = a + A a , the action of the air 
on the .airplane becomes F' , the resultant of F and AF, 



1 
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AF 1 being, if A a is very small, on the straight line & 
drawn through M parallel to the tangent at K to the polar p. 

From this it is deduced that the moment AM, which is pro- 
duced by the rotation AO< of the airplane, is equal to AF x c, 
c being the distance of the C.G. from the line d. 
We then have, at the limit for Aa= 0: 

AF ' 



and 



or, by approximation 



°s - zk x ° 



H = £ £r° 

TT A CO 



H 1 M X C 

F Aa 



and, if the resultant of K n and K y , or of K t and K n , is indi- 
cated by K: 



last ly 



K Aa 
h - 1 P x o 



or 



AK 

M 
jr 

3 



Aa „. 
h = lf - X c 



AE 

from which it follows that h is proportional to c, ^ being 
nearly constant for all angles practically accessible in flight 
and ~ being constant for each airplane. After this is estab- 
lished, it is possible to substitute, in the study of stability, 
the consideration of c for that of h, c having on its part 
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a more evident geometrioal and mechanical signif ication. To this 
o we have given the name "index stability. " 

We will see immediately that .the position of the line d, 
which depends on the angle of attack and the slope of flight, 
does not depend on the angle of deviation of the elevators. In 
fact, a variation At, of the angle of deviation f at constant 
speed, produces a variation of only AF p of the force of the air 
on the horizontal tail planes, which remains the same with refer- 
ence to the airplane, when the latter is tilted by the small an- 
gle Atx already considered. Consequently , both forces F and F« 
are displaced by combining with the same force AF p , but their 
geometrical difference does not vary. 

Consequently, not only the lines d can all be determined 
by a single position of the elevators, no matter whether the lat- 
ter is compatible or not with the equilibrium at the different 
angles of attack, but they remain unchanged for every slope of 
flight and angle of attack, even if the position of the C.G. or 
the axis of traction is changed. In fact, the displacement of the 
C.G. or the axis of traction has no other effect on the whole air- 
foil of the airplane than a variation of the angle of deviation 
of the elevators, which has no influence on the lines d. 

The lines d (each of which corresponds to a given value of 
r cc and of the slight slope <p and the distance of which from the 
C.G. is the index of stability) are therefore very useful in the 
study of the stability of an airplane, even before knowing the 
position of the C. G. , and can be utilized in this case for deter- 
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mining the most suitable location fbi* this point* 

* 

The construction of all the lines d is very simple in ac- 
cordance with the results of the tunnel tests of the complete air- 
plane or of just the cell and tail. In fact, the direction of 
each straight line is that of the tangent to the relative polar 
through the point which corresponds to the angle of attack under 
consideration. The distance of the intersection of the line d 
with the axis of the airplane from the point of this axis, with 
reference to which the moments are measured, is given by the ex- 
pression 

in which M is the coefficient of moment with reference to the 
^chosen point, K n is the component following the normal to the 
axis of the coefficient of total thrust K; and I is the wing 
chord. 

Now, M is composed of a part M c . due to the action of 
the air on the cell, and a part Mp, due to the tail planes. 
We may write 

S 

M P = xK np g£ t 

in which K np is the coefficient of normal thrust relative to the 

horizontal tail planes (function of * and a); S p the surface 

of the tail planes; t the distance (which may be considered 

constant) of the thrust K n p from the reference point of the mo- 
*We have given to the lines d the name of "metacentric lines. " 
In fact, each one is the locus of the metacenters of the airplane 
corresponding to each slope and angle of attack for different 
values of t. 
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ments. x is the ratio of the square of the velocity of the air 
with relation to the tail planes at the speed V of the airplane, 
x only needs to be considered when the tail planes are in the 

slip stream from the propeller and it is a determined function of 

V- P 

-J- or perhaps of ^g-^, P being the propeller thrust , or again 
of Kx + Ky tan <p, i. e. , an increasing function of the angle of 
attack and the slope of flight. We may write, according to 
Froude's theorem, 

S K x + K v tan® 
X ss 1 + — 



2 

in which S' e is a suitable fraction of the surface of the pro- 
peller and a 0 the density of the air. We have 

6K ^ 
x e 5£_ 

c' = I + ^ t 

dK n dK^ S 

d.,a 

6 K 

in which is dependent on the shape of the tail planes, a 

p 

constant which is equal to the complement to the unity of the 
ratio, nearly constant, of the angle: of deviation of the air cur- 
rent of the cell to the angle of attack a 0 f the cell, measured 
from the angle corresponding to the zero thrust. 

The result is that the sum of a term c'„ = -f-i which de- 

d K n 

pends only on the cell and of a term 



n € X~rr S 

q I „ 6 3 L. —E t 
P " £Kn_ S 

d-ot 

which, by x constant, is constant at all angles of attack and de- 
pends on the relation of certain chas'e.ct eristics of the tail 
planes to the corresponding characteristics of the cell, so that 
all the lines d relating to the airplane: are obtained on dis- 
placing, by a determined amount og.. (conatant, if x is cone-.:, 
st ant) , all the lines d relative to the single cell. 

Nothing is therefore more simple than the construction of 
the lines d, especially if we are not limited first to the con- 
sideration of the values of x equal to 1, i.e. , if we do not 
neglect the increase in efficacy of the tail planes due to the 
slip stream, which leaves us to consider first only the gliding 
stability. 

When we have the diagram of the lines d, nothing is sir^le 

4 

than the study of the longitudinal stability of the airplane. It 
is only necessary to measure the distance of the C. G. from the 
lines d corresponding to the conditions of flight, in order to 
obtain the values of the "index of stability. " 

If we cannot place the C. G. in a suitable position and if we 
are obliged to displace the lines d, the established formulas 
give very simply the variations of the characteristics of the 
tail planes neoessary for obtaining this displacement. 

Figs. 3 and 3 show the lines d for x = 1 for a separate 
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wing (Fig. 3) and for a complete airplane (Fig. 3)« 

This being established, there also, results a notable simpli- 
fication in the interpretation of the tests which we proposed in 
1919 for the verification, in full flight, of the calculations 
of stability, if we apply the latter to the determination of the 
4uantities c. 

In fact, if we consider two epeeds very near stable flight 
determined by the angles of attack a and a + Aa and by the 
slopes qp and cp + A cp, we have in the second speed the follow- 
ing variations of the forces (with reference to the unity of 
speed and the surface of the airplanes): the force AK acting ... 
along the line d, the force Ag^g = A ( K x + % tan cp ) acting 

on the axis of the propeller, at the distance e from the G.G. , 

S p 

the force cj K np Ax acting on the tail planes. 

If, as in the general case, these forces do not offset each 

other, it is necessary in the second case to deflect the eleva- 

6 K S 

tors by an angle at, so as to produce a force £ ^ p -g 35 cit 
whose moment will offset that of the forces considered above, i.e. 
it is necessary that 

6E n p Sp g 
- t x -gr At „ g A.-K + e A (K x + K y tancp) + t^B 

whence, at the limit, and on putting AK m AK y 
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If, therefore, we observe in flights at certain speeds, the 
values of t and if we determine the corresponding values of 
Ka£ + Kf tancp, K n p and x, by choosing a series of flights which 
may be considered as continuous, e. g. , the series corresponding 
to constant values of (Fig. 4) , we can calculate, by very sim- 
ple processes of graphic derivation, the values of o correspond- 
ing to each flight speed. 

If we choose the series for which a is constant, we have 
AK = 0 and 

A (K x + K y tan <p) - K y Aq> 

- t x |E a t e K v Acp + t |E K np Ax 

6 t S y S p 

whence, at the limit, 

*s Kn p = dS v tx s d9 e K y) 

which renders it possible by analogous methods to evaluate the 
term t Knp which would be more difficult by direct methods. 
If, according to Froude' s theorem, which is admissible, we 

put 

S K x + K y tanc P 
x ~ ^ + SP*- " a o 

• e T 

we have 

A(K Y + K, r tan d) 
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which enables the following very considerable simplification of 
the formulas established above 

c _ x ±hm fa t *L . / e e + East _ s p t ^ d(K x +K y ^jgj 

6t 8 dJGy A * -& 0 S ' e / "dK 

and for a constant 

lop ±b ^ Sg np d t 

2 Ky 

After many and long attempts, which have taken considerable U 
time and energy for perfecting the instruments and the general 
plan of the experiments, we have begun in Italy the study and sys- 
tematic determination of the characteristics of longitudinal sta- 
bility of airplanes in full flight, according to the method out- 
lined above. 

It is very probable that from these experiments, which will 
be tried also with airplanes known not to be very satisfactory in 
piloting, we will be able to draw Quite clear conclusions concern- 
ing the values whioh are, for the different types, the most favor- 
able for the good behavior of airplanes in flight and on the com- 
parison of values obtained by tunnel tests. 

Calculations of longitudinal stability of airplanes can then 
be made with confidence and profit, because we shall have eliminat- 
ed and clarified everything that has hitherto been complicated, 
obscure and indeterminate. 

Translated by the National Advisory Committee for Aeronautics. 
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